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Abstract 

In this paper we derive an expression for the static electric polarizability of a particle bound by 
a finite potential well without the explicit use of the continuum states in our calculations. This 
will be accomplished by employing the elegant Dalgarno-Lewis perturbative technique. 
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I. INTRODUCTION 



The study of the finite potential well model and its applications has played a significant 
part in our fundamental education for a long time. The model provides us with the educa- 
tional experience of studying a spectrum of wavefunctions which includes both bound and 
continuum (unbound) states. Second, by taking the appropriate limits, the finite potential 
model connects with the models of the infinite potential well and the attractive Dirac-delta 
potential. This approach helps us in relating the features of the different models, and pro- 
vides us with a useful approach for checking our results. Third, the bound states of the well 
produce the probability of the existence of the particle in the classically forbidden region 
(region where the total energy of the system is less than its potential energy). This feature 
helps with clarifying the differences between classical mechanics and quantum mechanics and 
we will study in our work the significance of the wavefunction in the classically forbidden 
region in determining the electric polarizability. 

The applications of the finite potential well are numerous and exist in many branches of 
physics. A well known example of a basic application is the modeling of the nucleon-nucleon 
interaction by a finite potential well. Another example is the modeling of the scattering of 
low energy electrons from atoms by the scattering of a beam of particles by a system bound 
by a potential well (Ramsauer Effect). - 

In our work we provide an example of the effectiveness and the simplicity of the method 
of Dalgarno and Lewis in determining the energy shift of a quantum level in second order 
perturbation theory- The energy shift will be used to calculate the electric polarizability 
due to the interaction between a static electric field and a charged particle which is bound 
by a finite potential well. As is known, the conventional method used to calculate the energy 
shift in the second order involves an infinite sum or an integral that contains all possible 
states allowed by the transition. 3 Some of these states, for example, scattering (unbound) 
states, can be very difficult or impossible to obtain in a large number of problems. The 
unperturbed state will be the only part of the system which we need to know to calculate 
the exact energy of that particular state to the second order when we apply the technique 
introduced by the Dalgarno and Lewis method. Following our approach helps in avoiding 
unnecessary approximations and mathematical difficulties. 3,4 This in turn allows us to pay 
more attention to the physics of the problem and to make our work beneficial to both 
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advanced undergraduate and first-year graduate students. Second, some of the results which 
we will obtain for the electric polarizability can not be reached if the conventional method 
is used. Third, we provide one more illustration of how simple methods can aid in our 
learning of perturbation theory. Fourth, in addition to our previous work on both statio^ 
and dynamic^ electric polarizability, we are providing one more example to demonstrate 
that the knowledge of the continuum is not always required to calculate the energy shift of 
a quantum level to the second order. 

In the next section, we give a summary of the method of Dalgarno and Lewis. We 
follow by introducing the necessary models, deriving the expression for the static electric 
polarizability and analyzing our results. In Section IIV[ we give some concluding remarks. 

II. METHOD OF DALGARNO AND LEWIS 

The method of Dalgarno and Lewis is based on replacing the conventional algorithm 
used to calculate the energy shift in second order perturbation theory with the solution of 
an inhomogeneous differential equation. We give in this section a summary of the method. 
References 2 and 7, Schiff 8 and Merzbacher- can be checked for more details. The ground 
state energy shift to the second order AE is given by 

_^ (iP \H'M^ n \H>\^ ) 

n=l CJ ° 

ip is the ground state and is occupied by a charged particle as we assume in our problem. E Q 
is the ground state energy. The functions ip n 's represent all states allowed by the transition 
due to the interaction H'. The summation in Eq. ([1]) excludes the state ip . 

The first step in eliminating the infinite sum (Eq. ([I])) is to find an operator F which 
satisfies the following relationship: 

%!^ = <*.|F|*>. (2) 

With the use of the completeness of all the states ^ n 's and simple manipulations, AE Q 
can be written as 

AE = (^ |F'F|0 - (ip \H'\il; )(iJj \F\ip o ). (3) 
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Now, with the use of the following property: 



(^1 [F, H ] |</> > = (E - E n )(i; n \F\^ ) 



(4) 



Eq. 02]), can be written as- 



(E - H )(j> = H'% - (^IF'IVo)^ 



(5) 



where 



H ^ = E 4>, 
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(6) 



and 



= F^, 



(7) 



It should be noted that <f> (Eq. (jSJ) can contain an arbitrary function x with % satisfying 
the relationship (i? — H )x = 0. The final expression of is taken to be orthogonal to 
^o- 1 ^- and with this choice Eq. ([3]) can be written as 



To obtain AE a using the method of Dalgarno and Lewis, we first find <p with Eq. (jSJ) 
being our starting point. The second and final step is to use in Eq. ([8]) to obtain AE Q . 
With these two steps, the only stationary state of the system we use in our calculation is ip Q 
and we completely avoid the infinite summation in Eq. (CQ). 

Before closing this section, we point to an important application of the Dalgarno and 
Lewis method. In our previous work of solving for the electric polarizability in the case of 
the delta potential,- 1 ^ \ nas been equal to zero. In the current work, and as we will find out, 
X has to exist to obtain the correct polarizability. 

III. THE ELECTRIC POLARIZABILITY 

As stated before, the model we are presenting in this article is of a particle bound by a 
finite potential well. The potential energy V(x) is given by: 



AE =^ \H'\ct>). 



(8) 
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-V (-a<x<a), 
V(x) = <( (9) 

(otherwise), 

where 2a is the width of the well. The one- dimensional Schrodinger equation with the given 
V(x) produces a spectrum of bound states in two classes: even parity states and odd parity 
states in addition to a continuum of unbound states.— The unperturbed state ip in our 
problem is the lowest energy even parity state for a given V Q and a. ip is then given by^ 



ip (x <-a)=N cos K ae ko{x+a) , (lO.a) 

and 

ipo(— a < x < a) = N cosK Q x, (10.6) 

and 

ifj Q (x >a) = NcosK ae- k ° {x - a \ (10. c) 

where E Q is the eigenenergy of the state ip a and is negative, K 2 = j£(V — \E \), k 2 Q = 2m r |f° 
and K Q a < |. The normalization constant N is given by 

N=- l - 7T- (11) 

I sin K a cos K a < cos 2 K a a 2 

The determination of the eigenvalues of the well bound states is conventionally obtained via 
finding the roots of certain transcendental equations.— For a given V and a, the transcen- 
dental equations needed to determine E Q are given by 

7 tan7 = /3 , (12) 

and 

ll + Pl = R\ (13) 

where 7 D = K Q a, (3 Q = k Q a and R 2 = 2m f 2 Vo . 

The conventional method as it appears in many textbooks,-^— does not provide the 
simplest procedure for determining the bound states' eigenvalues. Due to this a variety of 
useful approximations had been developed to replace the conventional method.— _ — In Ref. 



13, the authors used their approach of simplifying Eqs. (|T2|) and (|T3|) to approximate the 
finite potential well by a wider infinite potential well. The width of the infinite potential well 
is equal to the width of the finite potential well multiplied by the factor (1 + -^). The resulting 
approximate energy eigenvalues range in their accuracy with the ground state energy being 
the most accurate (less than 1% error) in the range of R close to or larger than four.— The 
results of Ref. 13, which we will get back to, will be useful in our work when we test our 
expression of the electric polarizability in the limiting case of large R. 

To determine the energy shift AE Q , we study the interaction between a particle of charge 
q that is initially in the state ip a and an external static electric field e which is applied in the 
vicinity of the particle. Once we find AE a , the electric polarizability a can be determined 
from the relationship AE Q = — ^ae 2 . Now, to simplify the geometry of the problem without 
any loss of generality, we take the electric field e to be parallel to the x axis. The electric 
dipole Hamiltonian which represents the interaction of the electric field and the charge is 
then given by 

H' = -qex. (14) 

At this point we are in a position to calculate AE Q and the polarizability. With H Q = 
11m" + V{x) where m is the mass of the particle and with the aid of Eqs. (10. a) and 
(10. c), we can write Eq. (jSJ) for the region x > \a\ as follows 

(f) ( x < -a) = -qexN cos K ae k ° {x+a) , (15) 

(j) ( x > a) = -qexN cos K ae- k ° {x - a) , (16) 

where <fi (x < —a) and (p(x > a) refer to the expressions of in the regions x < —a and 
x > a respectively. The second term on the right-hand side of Eq. (jSJ) vanishes because H' 
has odd parity and tfj has a definite parity. 

Equations (|15p and ffTB 7 ]) are very simple and can be solved essentially by inspection. The 
expressions for (x < —a) and <p(x > a) are given by 

N cos K ae k ° (x+a) , (17) 



-h 2 kl h 2 
+ 



2m 



d 2 



2m V dx 2 



and 



(J 2 



-h 2 k 2 h 2 

H 

2m 2m \ dx 2 



[x < 



( mqe 

\2h 2 k 



( mqe 



x2 + 1 ' 
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and 



(x > a) 



f mqe 

\2h 2 k 



( mqe 



N cos K ae 



-k (x—a) 



(18) 



Substituting the expressions of <p (x > a), ifi (x > a), <fi (x < —a) and ifi (x < —a) in Eq. (JHJ), 
and performing the necessary integration, we obtain the energy shift corresponding to if) in 
the region x > \a\. The electric polarizability, «i, corresponding to this value of the energy 
shift can be written as 



Oi 



mq 2 N^ 

h 2 



cos 2 K Q a 



a 5a 5a 5 
l$ + 2kJ + W + 4k E OJ 



(19) 



By taking the infinite potential well limit {K Q a — > ~) , a\ becomes zero. Now by taking the 
limits of 2a approaching 0, with V Q approaching infinity while keeping aV a constant, we can 
obtain a\ for the case of the attractive delta. For the given limits, iV 2 cos 2 K a will become 



k Q and «i is then given by | 
case of the attractive delta." 



[mq ) 
WW) 



which is the correct result for the polarizability in the 



+ V 



-a < x < a) 



-qexN cos K D x, 



(20) 



To evaluate the contribution to AE Q from the region —a < x < a, we first solve Eq. §5§ 
for (p in this region. Eq. §5§ for this region is given by 

' h 2 k 2 \ t^&_ 

2m J 2m dx 2 

where t is a trial function which will satisfy Eq. f[2TJ|) . <p t can be found by a simple process 
of trial and error and it is given by 

-mqeN\ 2 ._ T ^ f —mqsN\ 



-a < x < a) 



x sin K x + 



2h 2 K 2 J 



x cos K x. 



(21) 



2h 2 K Q J 

The expressions of (f> t (—a < x < a) and if> (—a < x < a) are then used in Eq. flS} to find the 
energy shift corresponding to if) in the region —a < x < a. The electric polarizability (a^)* 
corresponding to this energy shift is given by 



mq 2 N^ 

h 2 



— a^ a? 

^H + ^^s2K a 



5a 2 

— sm2K a 



cos 2K Q a H — sin 2K Q a 

8Kx 



(22) 



By taking the limit corresponding to the attractive delta, {pt?) t becomes zero. By taking 
the infinite potential well limit (K Q a -> f ), (a 2 ) f = -0.1324176 (^-^j- At this point the 
electric polarizability a is equal to \a.\ + (0:2) J. «i gives zero for the case of the infinite 
potential well as expected since ifj for the region x > \a\ represents the system outside 
the well. The contribution to a for the case of the infinite potential well then comes from 



(ot2) t an< i this means that a is negative. Checking the expression of AE Q (Eq. and the 
definition of a (AE Q = — |ae 2 ), we conclude that a should be positive. So, even though the 
expression of <p t (—a < x < a) satisfies Eq. ( 1201) . it can not be a complete expression for <p in 
the region —a < x < a. 

As we learned from the previous section, <fi can contain an arbitrary function x with the 
condition (E a — H D ) x = 0. We take x to be equal to ( -^r^ J C sin _K" D a; with the coefficient 
C being x independent. The reasons for making this choice for x are: (1) (0< + x) is expected 
to be orthogonal to ip -^^-'— (2) x should provide positive contribution to a; (3) The structure 
of x should be consistent with the structure of (f> t - 

To determine the coefficient C, we follow two simple steps. The first step is to calcu- 
late a for the infinite potential well with the use of the Dalgar no- Lewis technique and the 
conventional method. The comparison of the two results will then produce Ci n f, po t. (C in 
the case of the infinite well). In applying the Dalgarno- Lewis technique we calculate a by 
replacing (j) by <f) t + x in Eq. ([8]), setting K Q = an d performing the necessary integration. 
In using the conventional method, we have to use Eq. (CQ), however the equation contains 
an infinite sum. The infinite sum causes no problem, since the numerator inside the sum 
decreases rapidly with increasing n, and the denominator is proportional to (n 2 — 1). A sim- 
ple calculation using the spectrum of the infinite well will show that the transition n = 1 to 
n = 2 is far superior in magnitude to the next few terms. Using only the transition n — 1 to 
n = 2, and comparing the results of both the conventional and the Dalgarno-Lewis method, 
we find that the choice Ci n f, pot . = —a 2 produces an almost perfect agreement. The one- 
term conventional method produces a polarizability of 0.0701371 ^ m ^ 2 a ^ while substituting 
by Qn/.pot. = —a 2 i n the other expression of the polarizability (Dalgarno-Lewis) produces 
otinf.pat. which is given by: 

(2 4 \ 
^|^J . (23) 

The second step is to use the result Ci n f, pot . to help us predict C. Our prediction should 

produce the right result for the infinite potential well, satisfy the limiting case of the attrac- 

~(-) 2 

tive delta potential and be consistent with the structure of 4>f The choice of C = V-V 
clearly fulfills the first and third conditions and its fulfillment of the second condition will 
be demonstrated. In addition testing our choice of C away from the limit K a — > | is 
warranted. 
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At this point and with our choice of C, <j) ( — a < x < a) and the corresponding electric 
polarizability a 2 are given by 

, 2 



-a < x < a) 



-mqeN\ 
2h 2 K ) 



x 



x sin K Q x + — cos K Q x 



(-) 

sin if z 



and 



where 



and 



a 2 



mq 2 N^ 

h 2 



—a 
3K 2 



+ fi (K Q , a) cos 2K Q a + f 2 (K a , a) sin 2K Q a 



fi(K o ,0) 



f2{K o ,0) 



2K 2 



5a 



-5a 2 5 



2Kt ' 
(f) 2 



(24) 



(25) 



(26) 



(27) 



4^3 

By taking the limits corresponding to the attractive delta, a 2 becomes zero. We take the 
electric polarizability a to be equal to ol\ + a 2 . 

To test our choice of C away from the infinite well limit, we use the results of Ref. (13) 
to write an approximate expression for the polarizability (a apr ) and use this expression to 
compare with our expression of a. According to the work of Ref. (13), a finite potential 
well of width 2a can be approximated by a wider infinite potential well of with 2b with b = 
(l + -=;) a. The approximation produces small error in the eigenvalues and the probability 
densities of the eigenfunctions belonging to the lowest energy states for R close to or larger 
than four.— For example, in comparing the approximated and the exact energy eigenvalue 
for the lowest two energy states at R — 4, it is 0.69% for the lowest energy even parity 
state, and 3.15% for the lowest energy odd parity state. The approximation also produces 
very small error for the probability densities of the two states as a function of x (Fig. 3 
of Ref. (13)). These results are significant to our intended comparison due to two main 
reasons. First, the lowest energy state is the only state of the system which we use in our 
calculation of the polarizability. Second, the transition n = 1 to n = 2 produces almost all 
the contribution to the electric polarizability in the case of the infinite potential well. 



To present our comparison, we first write a apr by replacing a in Eq. ( |23l) by b. a 



l apr 



IS 



then given by 



a 



apr 



0.0702247 



1 



(28) 
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Now we set g = {^jp-^J a 4 and divide a±, a<i and a apr by g to get the dimensionless polariz- 
abilities a[, a' 2 and a' apr . a' x and a' 2 will be given by 



and 



where 



and 



Q, 



N' 2 cos z 7o 



' 1 



5 5 5 



ft 2ft 2ft 4y9|J 



iV 



/2 



T~2 + ^ 7 °) cos 27 ° + & (To) sin 2 7o 
37, 



ff ^ 1 5 
7l J ~*i8 4 7 4 



(il! 

27 4 
(f) 2 



—5 5 



Accordingly iV /2 is written as 



(29) 



(30) 



(31) 



(32) 



JY' 



sin 7q cos 7 , cos 2 7 P 



(33) 



and a' = a[ + a' 2 - In Table [Q we present the numerical values for a[, a' 2 , ol and a' apr 
for R values ranging from around four to R becoming very large. This is the range of R 
where the results for a apr are reliable according to the previous discussion. Also, in this 
range of R, a' 2 is close or approximately equal to a' and this provides a crucial testing to 
our choice of C. Now in comparing the numerical results for a' and a' apr we find that we 
have an excellent agreement. Of course, we did not expect a perfect agreement because 



a 



apr 



is based on an excellent approximation rather than exact calculations. To see how our 



choice of C contributes to the agreement between the numerical results of a' and a' apr , we 



define T as 



with (at, 



2Jt 



( Q a)t 
9 



T can be obtained by calculating the term 



■ cos 2 7o + sin 2 7o 



2 ^^^wo^ /0 i ^5-ojj._ „ and dividing the outcome by the bracket on the R.H.S. of 
q. ( 1301) for a given 7o . For 7o = 0.397T, we obtain T = 2.52, and with 7o = 0.47%, T = 2.84. 
This demonstrates that the choice of C is crucial to the obtained value of a' 2 (02). With this 
result, the observed dominance of a' 2 over a[ in determining a' (Tabled]) and the agreement 
of the numerical result of a' and a', we have verified the validity of our choice of C. The 
result of a = «i + «2 is then our final expression for the electric polarizability of the finite 
well. 
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Referring to Table HT1 we present our results for a[, a' 2 and a' in the region of 7 G < 0.2tt. 
In this region a[ which has been obtained with the aid of ip (x > \a\), is approximately equal 
to a'. In addition, we already verified that the delta-potential polarizability is obtained from 
ai when the appropriate limits are taken. These results demonstrate the superiority of the 
bound state in the classically forbidden region in determining the electric polarizability in 
the region of smaller j . 

Considering the dependence of a' on R, our results in both tables in addition to the 
results which are not shown here demonstrate that as R decreases, a! increases. With our 
definition of a' (a' — j), we fixed a and this means that a smaller R leads to a lower binding 
energy, i.e. lower \E D \. Since the polarizability is inversely proportional to the energy gap 
between the ground state and the higher energy states (Eq. (TjQ) and the basic definition of 
polarizability), we then expect a' to be increasing while R is decreasing. 

IV. CONCLUSION 

In this paper we have derived an expression for the static electric polarizability of a 
particle bound by a finite potential well. The only stationary state of the system which we 
have used in our calculation is the unperturbed state of the particle. This result and the 
studies appearing in our previous work,- 1 ^ teach us that the expression of the continuum 
is not required for every calculation of the electric polarizability. It is the elegance of the 
Dalgarno-Lewis method which has allowed us to avoid any complications arising from dealing 
with the continuum and to have only simple functions (sines and cosines) appearing in our 
expressions. The simplicity of the method and expressions helps the reader in checking every 
step of our work. It aids in clarifying the physics of the problem and certainly encourages 
physics students to learn about perturbation theory via the application of the Dalgarno- 
Lewis method. 

In deriving the expression for the static electric polarizability in the case of the delta- 
function,-^ the solution of Eq. (J5J) gave us the final expression for (j). In the case of the finite 
potential well a function \ has to be added to the solution of Eq. (J5}. In determining we 
used basic physics principles, the model of the infinite potential well, and the approximation 
of Ref (13). This process in turn has provided us with some valuable lessons. First, the 
use of the Dalgarno-Lewis method in calculating the electric polarizability can be extended 
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beyond problems which commonly appear in literature. Among such problems are the cases 
of the single delta, the infinite potential well, and the hydrogen atom in its ground state. 
Second, we avoided the infinite sum of the conventional method, but we still used the first 
term of the sum in the case of the infinite potential well to determine x- 111 this process, we 
continue to depend only on bound states in the simplest form to obtain our result. Third, 
it is the realization that the approximation of Ref. (13) has practical applications in solving 
problems related to the finite potential well. 

The starting point in our derivation has been the solution of Eq. (jSJ). This allowed us 
to write the contribution of a in terms of ai and ol-i- This separation has been helpful in 
clarifying how the bound state in the classically forbidden region contributes to the electric 
polarizability. In addition, we can specify that a,\ is what produces the limit for the delta 
potential while «2 produces the limit for the infinite potential well. Such separation would 
be impossible in the case of the conventional method since the determination of a depends 
on the calculation of the term |(-0„|if' |-0 O )| 2 (Eq. (TTJ). 
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7o 


Po 


R 


«i 


a' 2 


a' 


a' 

^apr 


0.39vr 


3.403183 


3.617018 


0.015178 


0.173148 


0.188326 


0.186438 


0.417T 


4.433507 


4.616825 


0.005510 


0.147482 


0.152993 


0.153844 


0.43vr 


6.043511 


6.192650 


0.001663 


0.125180 


0.126843 


0.127803 


0.45vr 


8.925856 


9.037118 


0.000363 


0.106019 


0.106382 


0.106858 


0.47vr 


15.620252 


15.589884 


3.99E-5 


0.089754 


0.089794 


0.089913 


0.49vr 


48.983879 


49.008061 


4.24E-7 


0.076129 


0.076129 


0.076134 



TABLE I: Comparison of the dimensionless polarizability a' (our calculation with the Dalgarno- 
Lewis method) with the dimensionless polarizability a' apr (our calculation with the infinite potential 
well approximation of Ref. (13)). 



7o j3 R a[ a' 2 a' 

0.19vr 0.405655 0.721698 4.93E+1 0.620993 4.99E+1 

0.17vr 0.315849 0.620477 1.31E+2 0.677762 1.32E+2 

0.15vr 0.240108 0.528884 3.87E+2 0.733438 3.88E+2 

TABLE II: Comparison of the contributions of a[ and a! 2 to a' for j < 0.2ir. 
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